Magnetic Properties of J-J-J' Quantum Heisenberg Chains with Spin S = 1/2, 1, 3/2 
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By means of the density matrix renormalization group (DMRG) method, the magnetic properties 
of the J-J-J' quantum Heisenberg chains with spin S = 1/2, 1, 3/2 and 2 in the ground states are 
investigated in the presence of a magnetic field. Two different cases are considered: (a) when J 
is antiferromagnetic and J' is ferromagnetic (i.e. the AF-AF-F chain), the system is a ferrimag- 
net. The plateaus of the magnetization are observed. It is found that the width of the plateaus 
decreases with increasing the ferromagnetic coupling, and disappears when J'/J passes over a crit- 
ical value. The saturated field is observed to be independent of the ferromagnetic coupling; (b) 
when J is ferromagnetic and J' is antiferromagnetic (i.e. the F-F-AF chain), the system becomes 
an antiferromagnet. The plateaus of the magnetization are also seen. The width of the plateaus 
decreases with decreasing the antiferromagnetic coupling, and disappears when J' /J passes over a 
critical value. Though the ground state properties are quite different, the magnetization plateaus 
in both cases tend to disappear when the ferromagnetic coupling becomes more dominant. Besides, 
no fundamental difference between the systems with spin half-integer and integer has been found. 
PACS number(s): 75.10.Jm, 75.30.Kz 
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I. INTRODUCTION 

Low-dimensional quantum spin systems have been at- 
tracting both experimental and theoretical interest in the 
last decades due to an interplay of strong quantum fluc- 
tuations and topology. Several theoretical predictions for 
the low-dimensional quantum spin chains have been veri- 
fied by experimental studies. For the Heisenberg antifer- 
romagnetic (HAF) chains with spin S — half — integer, 
the celebrated Lieb, Schultz and Mattis theorem showed 
that the excitation of the system is gapless 0. For 
the HAF chains with spin S — integer, the excita- 
tion from the singlet ground state to the triplet excited 
state was conjectured to be gapful, now known as Hal- 
dane conjecture 3]. Another interesting phenomenon in 
HAF spin chains is the occurrence of the magnetization 
plateaus, that can be viewed as an essentially macro- 
scopic quantum phenomenon, and has gained much at- 
tention recently. A decade ago, Hida has considered an 
S = 1/2 HAF chain with exchange coupling of 3-site 
translational invariancc in the presence of an applied 
magnetic field, and uncovered a plateau in the magne- 
tization curve at 1/3 of the saturation magnetization 
Slightly after Hida's numerical calculation on the plateau 
of the S = 1/2 F-F-AF chain, an analytical approach was 
done by Okamoto 5j. These results lead to a more general 
necessary condition for the appearance of the magnetiza- 
tion plateaus which was proved by Oshikawa, Yamanaka 
and Affleck 6] . It tells us that for the HAF spin chains 
with S = integer or half — integer, the magnetization 
curve can have plateaus at which the magnetization per 
site m is topologically quantized by 



n(S — m) — integer, 



(1) 



where S is the magnitude of the spin, and n is the period 
of the ground state determined by the explicit spatial 
structure of Hamiltonian. At the plateaus, the spin gaps 
open, that can be in some sense regarded as a kind of 
generalization of the Haldane conjecture. Similar to the 
quantum Hall effect, the magnetization plateau is an- 
other striking example of the macroscopic quantum phe- 
nomenon, in which magnetization is quantized to frac- 
tional values of the saturated magnetization value and is 
a function of the magnetic field. 

The magnetization plateaus are predicted and 
observed in many low-dimensional spin systems. 
Among others, a magnetization plateau at half 
the saturation magnetization was observed in 
5 = 1 HAF bond-alternating chain compounds 
[Ni 2 (dpt) 2 (M-ox)(/i-N 3 )](PF 6 ) (abbreviated as 
NDOAP) and [Ni(333 - tet)(/x - N0 2 )]C10 4 (abbre- 
viated as NTENP), where the experimental result is 
in agreement with the numerical calculations 0, S 
the ferrimagnetic mixed spin chains such as the 
bimetallic chain MM'(pbaOH)(H 2 0)nH 2 and the 
organic compound Mn(hfas) 2Q (3R )2 show quantum 
magnetization plateaus 0, ITU Tl2l Il3l IT^ : the mag- 
netization plateaus are also found in p-mcrized chains 
and ladders E E E3. The effect of randomness 
on magnetization process also attracts much theoretical 
interest [Tgt |2C| . One-dimensional (ID) helical spin 
system Co(hfac) 2 NITPhOMe shows some interesting 
magnetic behaviors, where one of the unusual properties 
is that the magnetization shows plateaus at zero and 
■| of the saturation if a magnetic field is applied along 
the helical axis, but no plateaus if the field is applied 
in the plane perpendicular to that axis|2lj. Another 
intriguing topic in ID spin systems is concerned with the 
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spin-orbital model, such as Na2Ti2Sb20 and NaV205, 
where the orbital degree of freedom plays an imp ortant 
role in the appearance of plateaus [23, 1221 12J, [23 ■ 

Despite the existence of considerable theoretical in- 
sight, a satisfactory understanding of the experimen- 
tal situation is still sparse. Experimentally, polynuclear 
azido-bridged derivatives are a rich source of new mag- 
netic systems. Because of the extreme versatility of the 
coordination modes of the azido ligand, these magnetic 
systems may coordinate as ^1,3 — N 3 (end-to-end, EE), 
Hi t i — N3(end-on, EO), or even in more exotic modes as 
Mm,i - N 3 or ^1,1,3 - N 3 [H I2II23. The similar de- 
gree of stability of the EE or EO coordination modes 
often favors a variety of topologies or dimensionalities. 
The antiferromagnetic (AF) interaction and the ferro- 
magnetic (F) coupling are generally found in the EE and 
EO mode, respectively. The exotic topologies with al- 
ternating patterns EE/EE/EO have been reported [2^| . 
which offers exciting research prospects. 

In this paper, motivated by the exotic magnetic prop- 
erties of the period n — 3 quantum spin system such as 
EE/EE/EO, we shall investigate the ground state of a 
ID J- J- J' quantum spin chain with spin S — 1/2, 1, 3/2, 
and 2, respectively. It is found that in the magnetic pro- 
cess, magnetization plateaus are observed, in agreement 
with Eq.QJ. The width of the magnetization plateaus is 
found to depend on the ratio of J'/ J. The fundamen- 
tal difference of the properties between the systems with 
spin half-integer and integer is not observed. 

The rest of this paper is outlined as follows. In Sec. 
II, we introduce the model Hamiltonian for the ID J-J- 
J' Heisenberg spin chain. In Sec. Ill, we present our 
numerical results of the ground state of the system. A 
brief summary is given in Sec. IV. 



(a) AF-AF-F 

(b) F-F-AF 



FIG. 1: The spin arrangements of the J- J- J' chain, where f 
denotes spin up, J, denotes spin down, J and J' are exchange 
couplings, (a) AF-AF-F chain, J = Jaf, J' — Jf', (b) F-F-AF 
chain, J — Jf, J' = Jaf- 

the system becomes a uniform Heisenberg ferromagnetic 
chain. It has no magnetization plateaus before saturation 
in this case. 

When J ^ J', the system shows complex behaviors. 
As indicated in Fig. ^ there are two interesting cases: 
(a) J = Jaf, J' — Jf- It is a ferrimagnetic spin chain (i.e. 
the configuration is AF-AF-F). (b) J = J F , J' = Jaf- It 
is an antiferromagnetic system (i.e. the configuration is 
F-F-AF). In order to probe the fundamental difference of 
properties of the systems with spin half-odd integer and 
integer during the magnetizing process, we shall consider 
S = 1/2, 1, 3/2 and 2 for each case below. 



III. DMRG RESULTS 



II. MODEL 

Motivated by the experimental study of the magnetic 
coordinated co mpoundsHHHEI, let us investigate the 
ID J-J-J' Heisenberg quantum spin system. The Hamil- 
tonian of the system reads 

H = JS3j_2'S3j_i + JS3j-i-S3j + J'Ssj rS^j+i) — h 

3 3 

(2) 

where J and J' are exchange integrals with J, J' > 
denoting the AF coupling and J, J' < the F coupling, 
h is the external magnetic field, and we take gag = E 
The schematic spin arrangement of the system is shown 
in Fig. [2 

When J = J' = Jaf, the system becomes a uniform 
HAF spin chain. In this case, it is well-known that for 
S = half — integer, there is no magnetization plateau in 
the magnetic process before saturation; for S = integer, 
owing to the existence of the Haldane gap there appears 
a plateau at zero magnetization below the lower critical 
field where the spin gap closes. When J = J' = J F , 



The density matrix renormalization group 
(DMRG)[3l| was proposed more than ten years ago, 
which nowadays becomes a powerful numerical method 
invoked to study the ground state and low-lying states 
of low-dimensional lattice systems. In the following, the 
magnetic properties of the ground state of the spin chain 
with configurations (a) and (b) in Fig. ^ with open 
Y2 boundary conditions are investigated by the DMRG 
3 'method. In our calculations, we took the chain length 
L = 60, the number of states kept per block N = 60 for 
spin S =1/2 and 1, and N = 80 for spin S = 3/2 and 2. 
The truncation error is less than 10~ 3 in all cases. 



A. AF-AF-F CHAIN 

Let us first consider the AF-AF-F chain, as shown in 
Fig. nja). In this case, J = Jaf and J' = J F . Obvi- 
ously, it is a ferrimagnetic chain. When Jp / Jaf is much 
less than 1, the antiferromagnetic coupling is dominant; 
while Jf/Jaf is much greater than 1, the ferromagnetic 
coupling becomes dominant. Without loss of generality, 
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AF-AF-F 




FIG. 2: Magnetic curves of the AF-AF-F chain of S = 1/2, 1, 
3/2 and 2, respectively, where h is the applied magnetic field, 
m is the magnetization per site. J = Jap = 1, J' = Jf- 

we shall take the ratio Jf/Jaf m the range of [0.1, 10] in 
this subsection. 

Fig. |21 show the magnetization process of the AF-AF- 
F chain with spin S — 1/2, 1, 3/2 and 2, respectively. 
For S = 1/2, the two plateaus at m — 1/6 and m — 
1/2 (saturation plateau) are obtained, consistent with the 
necessary condition given by Eq.Q with n = 3. For S — 
1, the three plateaus at m = 1/3, 2/3 and 1 (saturation 
plateau) are seen. For S = 3/2, the four plateaus at to = 
1/2, 5/6, 7/6 and 3/2 (saturation plateau) are observed, 
and for S = 2, the five plateaus at m = 2/3, 1, 4/3, 
5/3 and 2 are obtained. The number of the plateaus is 
(2S +1). It is clear that the smaller the ratio Jf/Jaf 
is, the more obvious the plateaus are, that appears to be 
independent of the magnitude of spin. In other words, 
when the AF interaction in the system plays a dominant 
role, the magnetization plateau becomes wider; when the 
ferromagnetic interaction is dominant, the width of the 
plateau gets narrower, as indicated in Fig. [3 In between 
the plateaus, the magnetization per site increases with 
increasing the field for all cases with different spins. In 
this system, the saturated field does not change with the 
ratio Jf/Jaf- A further discussion on this point can be 
found in Sec. IV. For the present system with different 
spins S = 1/2, 1, 3/2 and 2, the qualitatively similar 
behaviors of the magnetization process are observed. 

The appearance of the magnetization plateaus can be 
understood from the spatial dependence of the averaged 
local magnetic moment, as discussed in Ref. 30]. Fig. 
EH presents the spatial variation of (Sj) in the ground 
state for the AF-AF-F chain with spin S — 3/2 and 2, 
respectively. For both cases, the ratio Jf/Jaf = 0.1 
is taken. Fig. E^a) shows the case with S = 3/2. 
At the external field h/JAF = 2, the expectation val- 
ues (Sj) versus site j follow the sequence such that 
{...,(1.2268,1.2268,-0.9535),...}, resulting in the mag- 
netization per site to = Ylj=i({S£j-2) + (^3j-i) + 
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FIG. 3: The spatial variation of {Sj) in the ground state for 
the finite AF-AF-F chain with lengthL = 60 at Jf/Jaf = 0.1. 
(a) For S = 3/2, the external field is taken as /i/Jaf = 2, 3, 
4, respectively, (b) For S — 2, the external field is taken as 
h/JAF = 2.5, 3.5, 4.5 and 5.5, respectively. 

(S$j))/3N = 1/2, corresponding to the plateau m = 1/2. 
When the field increases to h/Jxe — 3, the behav- 
ior of {SI) becomes {..., (1.2280, 1.2280, 0.0440), ...}, giv- 
ing rise to m = 5/6; at H/Jaf — 4, {Sj) versus j 
behaviors as {...,(1.3390,1.3390,0.8219),...}, leading to 
to = 7/6. Fig. OUb) shows the case with S = 2. Sim- 
ilarly, at h/J A p — 2.5, the expectation value {Sj) ver- 
sus site j behaviors as {..., (1.7178, 1.7241, -1.4448), ...}, 
resulting in m = 2/3; at H/Jaf — 3.5, (Sj) 
varies according to {...,(1.6478,1.6478,-0.2956),...}, 
giving the m = 1; at H/Jaf — 4.5, (Sj) be- 
haviors as {...,(1.7137,1.7127,0.5746),...}, leading to 
to = 4/3; and at h/JAF = 5.5, it becomes 
{..., (1.8390, 1.8390, 1.3220), ...}, giving rise to to = 5/3. 
These to values just correspond to the magnetization 
plateaus. In addition, the spin configuration for the pa- 
rameters above are stable against the small local field 
h\oc/ Jaf—0-1, which is set at j = 1 site and parallel to 
the external magnetic field h. 

Fig. ^illustrates the phase diagram in the h c — Jp / Jaf 
plane for the AF-AF-F chain with spin S = 1/2, 1, 3/2 
and 2, respectively. In the shaded regions, the magne- 
tization plateaus appear, implying that the excitations 
from the ground state are gapful. When Jf/Jaf in- 
creases beyond the critical value (Jf/^af)c at which the 
plateau vanishes, i.e. the ferromagnetic coupling becomes 
more dominant, the magnetization plateaus tend to dis- 
appear. The plateau-non-plateau transition could be of 
the Berezinskii-Kosterlitz-Thouless(BKT)-type, and the 
width of the plateau to = Ah c near the transition point 
behaves as 

Ah c = Dexp(-C/v/(J^/J AF ) c - 3 F /3 AF ), (3) 

where D,C are constants[32|. 

As the width of the plateau might be extremely small 
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FIG. 4: (Color online) Phase diagram in the h c — Jf/Jaf 
plane for the AF-AF-F chain with S = 1/2, 1, 3/2 and 2, re- 
spectively. In the shaded regions, the magnetization plateaus 
labeled with the quantized values of m appear. 
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FIG. 6: Magnetization process of the F-F-AF chain with 
S = 1/2, 1, 3/2 and 2, respectively, where h is the ap- 
plied magnetic field, and m is the magnetization per site. 
J = Jf = 1, J — Jaf ■ 
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FIG. 5: The plateau-non-plateau transition point (Jf/Jaf)c 
for the AF-AF-F chain with S = 1, 3/2 and 2, respectively, 
where D, C are constants in Eq. and A/i c (Jf/Jaf) are 
scaled by the value A/i c (0.1). 



near the transition point, a direct estimation of the crit- 
ical point from the numerical data may be somewhat in- 
accurate. However, Eq. Q can be used to estimate the 
plateau-non-plateau transition point (Jf/ Jaf)c from the 
raw numerical data, which yields (Jf/Jaf)c = 4.5 for 
S = 1 and m = 2/3, 2.2 for S = 3/2 and m = 7/6, and 
1.2 for S = 2 and m = 5/3, respectively, as shown in Fig. 
|SJ The larger the spin S is, the smaller the critical value 
of (Jf/Jaf)c is. It is also noticeable that for a given S, 
the critical values for different plateaus are slightly dif- 
ferent. No fundamental difference between the systems 
with spin integer and half-integer is observed. 



B. F-F-AF CHAIN 

Now we consider the case with J — Jf, J' = Jaf, 
namely, the system is a F-F-AF chain, as shown in Fig. 
^b). It is an antiferromagnet. In contrast to the AF- 
AF-F chain, when Jaf/Jf is much less than 1, the ferro- 
magnetic coupling is dominant; while Jaf/Jf is greater 
than 1, the antiferromagnetic interaction predominates. 
We shall take the ratio Jaf/Jf in the range of [0.1, 8] in 
this subsection. 

Fig. shows the magnetization process of the F-F-AF 
chain with different spins. For S = 1/2, we have ob- 
served the two plateaus at m = 1/6 and 1/2, consistent 
with the condition of Eq.Q with n = 3. For 5=1, 
the three plateaus at m = 1/3, 2/3 and 1 (saturation 
plateau) are observed. For S = 3/2, the four plateaus 
at m = 1/2, 5/6, 7/6 and 3/2 (saturation plateau) are 
seen, and for S = 2, the five plateaus at m = 2/3, 1, 4/3, 
5/3 and 2 are obtained. It is seen that the smaller the 
ratio Jaf/Jf is, the less obvious the plateaus (except for 
the saturation plateaus) are. In other words, when the 
AF coupling is predominant in this system, the plateau 
is more obvious, and the width of the plateau is much 
wider. In between the plateaus, the magnetization in- 
creases rapidly with the magnetic field. From Fig. one 
may find that, the plateaus at m = for S = 1, m = 1/6 
for S = 3/2, and m — and 1/3 for S — 2, which are 
also allowed by Eq.Q) with n — 3, are not observed. 
This is because Eq.Q is a necessary condition for the 
occurrence of the magnetization plateaus. Therefore, the 
number of the plateaus is still (2S* +1). As the present 
system is an antiferromagnet, the magnetic curves look 
different to those of the AF-AF-F chain. The satura- 
tion fields depend closely on the magnitude of the ratio 
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FIG. 7: The spatial variation of in the ground state 

for the finite F-F-AF chain with length L — 60 at Jaf/Jf = 
8. (a) For S = 3/2, the external field is h/J F = 5, fO, 20, 
respectively, (b) For S — 2, the external field is /i/Jf = 5, 
f0, 20 and 25, respectively. 
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FIG. 8: (Color online) Phase diagram in the h c — Jaf/Jf 
plane for the F-F-AF chain with S = 1/2, 1, 3/2 and 2, re- 
spectively. In the shaded regions, the magnetization plateaus 
labeled with the quantized values of m appear. 



Jaf/Jf, in contrast to the AF-AF-F chain. In addition, 
the similar behaviors for the system with different spins 
are obtained, and no fundamental difference of the prop- 
erties of the system with spin integer and half-integer is 
observed. 

Fig. [7| shows the behavior of (S?) against j for 
the F-F-AF chain with spin 5 = 3/2 and 2, respec- 
tively. For both cases, the ratio Jaf/Jf — 8 is set. 
Fig. E[a) shows the case with S = 3/2. For the ex- 
ternal field h/Jp = 5, the spatial variation of (Sj) 
is {...,(0.0366,0.0365,1.4269),...}, giving rise to m = 
1/2. As the field is h/Jp = 10, the behavior of 
(S?) against j becomes {...,(0.5071,0.5071,1.4857),...}, 
resulting in to = 5/6. For h/J-p — 20, the spin 
configuration is {...,(1.0018,1.0018,1.4964),...}, giving 
to = 7/6. Fig. EJb) shows the case with S — 
2. For h/Jp — 5, the expectation values (Sj) is 
{...,(0.0871,0.0863,1.8265),...}, leading to m = 2/3. 
When h/Jp = 10, the behavior of (Sj) versus j be- 
comes {..., (0.5181, 0.5183, 1.9636), ...}, giving m = 1; at 
h/Jp = 20, it becomes {...,(1.0060,1.0060,1.9882),...}, 
yielding to = 4/3; at h/Jp — 25, the spin configura- 
tion is {..., (1.5018, 1.5018, 1.9964), ...}, giving to = 5/3. 
The to values are just corresponding to the magnetiza- 
tion plateaus. In addition, the spin configuration for the 
parameters above are stable against the small local field 
hioc/ 'Jf=0.1 at j = 1 site, which is parallel to the exter- 
nal magnetic field h. 

Fig. |S] presents the phase diagram in the h c — Jaf/Jf 
plane for the F-F-AF chain with spin S = 1/2, 1, 3/2 
and 2, respectively. In the shaded regions, the magneti- 
zation plateaus occur. At the plateaus, the spin excita- 
tions from the ground state are gapful. When Jaf/Jf in- 
creases beyond the critical value (Jaf / Jp)c at which the 
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FIG. 9: The plateau-non-plateau transition point (Jaf/Jf)c 
for the F-F-AF chain with S = 1/2, 1, 3/2 and 2, respectively, 
where D, C are constants in Eq. and A/i c (Jaf/Jf) are 
scaled by the value A/i c (1.0). 

plateau occurs, i.e. the AF coupling becomes more dom- 
inant, the magnetization plateaus tend to appear. The 
plateau-non-plateau transition is also expected to be of 
the BKT-type, and the width of the plateau m — Ah c 
near the transition point behaves as 

Ah c = Dexp(-C/VWJF - (WJf) c ), (4) 

where D, C are two constants [32|. 

As the width of the plateau might be extremely small 
near the transition point, Eq. Q is used to estimate the 
plateau-non-plateau transition point (Jaf/^f)c from the 
raw numerical data. It gives {Jaf / Jf)c — 0.09 for S = 
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1/2 and m = 1/6, 0.13 for 5 = 1 and m = 2/3, 0.18 for 
5 = 3/2 and m = 7/6, and 0.24 for 5 = 2 and m = 5/3, 
respectively, as shown in Fig. OH It should be pointed out 
that the breakdown of the magnetization plateau of 5 = 
1/2 F-F-AF chain has been studied, with the estimation 
of the critical value (Jaf/Jf)c — 0.065[22. It can be 
seen that the larger the spin 5 is, the greater the value 
of (Jaf/Jf)c is. For a given 5, the value of (Jaf/Jf)c 
for different plateaus is slightly different. In addition, 
no fundamental difference between the F-F-AF systems 
with spin integer and half-integer is also observed. 



IV. SUMMARY AND DISCUSSION 

In this paper, by using the DMRG method we have nu- 
merically studied the magnetic properties in the ground 
states of the J-J-J' trimerized quantum Heisenberg chains 
with spin 5 — 1/2, 1, 3/2 and 2, respectively. The two 
different cases are considered: (a) J = Jaf and J' = Jp, 
i.e., the AF-AF-F trimerized ferrimagnetic chain. The 
magnetization plateaus are observed when the AF cou- 
pling is dominated. The positions of these plateaus are 
allowed by Eq.QJ with n = 3, and the number of the 
plateaus is to be (25 + 1). For a certain spin 5 and a 
plateau (m ^ 5), the width of magnetization plateaus 
decreases with increasing the F coupling, and becomes 
zero when the ratio Jf/Jaf passes a critical value, e.g. 
( J f /Jaf)c = 4.5 for S = 1 and m = 2/3, 2.2 for S = 3/2 
and m = 7/6, and 1.2 for 5 = 2 and m = 5/3, re- 
spectively. The larger the spin 5 is, the smaller the 
critical value of (Jp/Jaf)c is- F° r a given 5, the value 
of (Jf/</af)c for different plateaus is slightly different. 
Furthermore, the saturation field dose not change with 
the F coupling, (b) J = J F and J' — Jaf, i.e., the 
F-F-AF trimerized antiferromagnetic chain. The criti- 
cal value in the plateau-non-plateau transition is about 
(Jaf/Jf)c = 0.09 for 5 = 1/2 and m = 1/6, 0.13 for 
5 = 1 and m = 2/3, 0.18 for 5 = 3/2 and m = 7/6, and 
0.24 for 5 — 2 and m = 5/3, respectively. For a given 5, 
the value of (Jaf/Jf)c for different plateaus is slightly 
different. We would like to point out that though the 
ground state properties of the F-F-AF chain are quite 
different from those of the AF-AF-F chain, as the former 
is a antiferromagnet, while the latter is an ferrimagnet, 
the magnetization plateaus in both cases tend to dis- 
appear when the ferromagnetic coupling becomes more 
dominant. 

As the width of the plateau might be extremely small 
near the transition point, a direct estimation of the criti- 
cal point from the numerical data may be somewhat inac- 
curate, as mentioned above. However, near the transition 
point the width of the plateau is expected to behave as 
Dexp(-C/y/\J'/J - (J'/J) c |), where D, C are two con- 
stants. According to this property, ( J'/ J) c can be esti- 
mated from the raw numerical data. Such an estimation 
is consistent with the previous study, showing that our 
present calculations are reliable. 



The appearance of the magnetization plateaus can be 
understood from the spatial dependence of the averaged 
local magnetic moment. For the AF-AF-F chain with 
spin 5 = 3/2, as the ratio Jf/Jaf = 0.1 and the exter- 
nal field h/ Jaf = 2 are taken, the spatial variation of 
(Sj) versus site j in the ground state follow the sequence 
such that {...,(1.2268,1.2268,-0.9535),...}, resulting in 
the magnetization per site m — ^jLi(( l ^'|j-2) + ('^3j-i) + 
(Sg -))/3iV = 1/2, corresponding to the plateau m = 1/2. 
Similar cases are found when the external field h / Jaf = 3 
and 4 are set. For the AF-AF-F chain with spin 5 = 2, 
F-F-AF chain with spin 5 = 3/2 and 2, the similar spa- 
tial dependence of the averaged local magnetic moment 
are observed. It should be noted that the spin config- 
uration for the parameters above are stable against the 
small local field /ii oc /Jaf=0.1 for AF-AF-F chain and 
7i1oc/Jf=0.1 for F-F-AF chain, which is set at j = 1 site 
and parallel to the external magnetic field h. 

By the spin-wave method starting from the fully polar- 
ized state, the saturation magnetic field of the dimerized 
and quadrumerized AF chains can be obtained analyti- 
cally [I3 • Inspired by this method, we have deduced an 
expression for the trimerized chains. The energy differ- 
ence AE between the lowest energy with M z — (3iV — 1)5 
and that with M z = 3NS is given by the lowest eigen- 
value of the following matrix, 

— (J' + 1) J'e lk e- tk \ 

J'e- lk -(J' + l) e lk 
e zk e -tk _ 2 J 

where k is the wave number, J = 1, J' and S are the 
same notation as mentioned above. We have found that 
the energy difference AE is the lowest when k = tt. 
This gives ri se to the satura tion field H s — ~AE = 
(3 + 2J' + v /9-4J' + 4(J') 2 )5/2. When J' = 1, it be- 
comes a simple AF chain, H s = 4S. When J' < , it 
becomes an AF-AF-F chain, and H s is nearly indepen- 
dent of J' . These results are also in agreement with the 
result of the quadrumerized chain [T^. Most importantly, 
all these analytic results are in agreement with our nu- 
merical results. 

In trimerized Heisenberg chains, the total number of 
the plateaus is 25 + 1, which implies that the OYA con- 
dition is only a necessary condition in trimerized chains. 
This observation can be understood straightforward in 
the case of ferrimagnet, i.e. the AF-AF-F chain. As 
shown in Fig. ^a), the possible minimum magnetization 
value is m = 5/3, so the plateaus whose magnetization 
value is smaller than 5/3 cannot naturally appear. For 
example, for 5 = 3/2, the plateaus permitted by the 
OYA condition should be m = 1/6, 1/2, 5/6, 6/7 and 
3/2, respectively. As the plateau m = 1/6 is smaller 
than the minimum magnetization of 5/3 = 1/2, it does 
not appear. Then, the number of the emerging plateaus 
is 4 = 25 + 1. However, this argumentation seems not 
to apply to the case of antiferromagnet, i.e. the F-F-AF 
chain. In Fig. ^b), the possible minimum magnetiza- 
tion value is m — 0. The plateaus m = 1/6 of 5 = 3/2 
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would be expected to emerge during the magnetization, 
but they do not appear. The reason why the plateau 
m = 1/6 of S = 3/2 vanishes is under exploration. It is 
worthy of stressing that for a certain spin S the number 
of disappearance of plateaus is smaller than that of the 
emergence of plateaus. 

The obtained results of the two cases reveal that the 
degree of the inhomogeneity of the couplings measured 
by the ratio J'/ J, not just the period of the quantum spin 
chain indicated in Eq.JTJ, determines its ground-state 
properties. In addition, no fundamental difference of the 
properties of the system with spin integer and half-integer 
is observed. It should be noted that the magnetization 
plateaus are only related to the AF coupling. In cases (a) 
and (b), when the F coupling is much stronger than the 



AF one, there is no magnetization plateaus at all, though 
the ratio of the couplings can still be large. Clearly, these 
magnetization plateaus come from the quantum origin, 
that are closely related to the competition between the 
exchange couplings. The thermodynamic properties of 
this system are studied under the way. 
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